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ABSTRACT 

We compute robust lower limits on the spin temperature, Tg, of the 2 = 8.4 intergalac¬ 
tic medium (IGM), implied by the upper limits on the 21-cm power spectrum recently 
measured by PAPER-64. Unlike previous studies which used a single epoch of reion¬ 
ization (EoR) model, our approach samples a large parameter space of EoR models: 
the dominant uncertainty when estimating constraints on Tg. Allowing Tg to be a free 
parameter and marginalizing over EoR parameters in our Markov Chain Monte Carlo 
code 21CMMC, we infer Tg ^ 3 K (corresponding approximately to lcr) for a mean IGM 
neutral fraction of xhi A) 0.1. We further improve on these limits by folding-in addi¬ 
tional EoR constraints based on: (i) the dark fraction in QSO spectra, which implies 
a strict upper limit of x\h\z = 5.9] ^ 0.06 + 0.05 (lcr); and (ii) the electron scattering 
optical depth, r e = 0.066 ± 0.016 (lcr) measured by the Planck satellite. By restricting 
the allowed EoR models, these additional observations tighten the approximate lcr 
lower limits on the spin temperature to Tg ^ 6 K. Thus, even such preliminary 21-cm 
observations begin to rule out extreme scenarios such as ‘cold reionization’, implying 
at least some prior heating of the IGM. The analysis framework developed here can be 
applied to upcoming 21-cm observations, thereby providing unique insights into the 
sources which heated and subsequently reionized the very early Universe. 

Key words: galaxies: high-redshift - intergalactic medium - cosmology: theory - 
dark ages, reionization, first stars - diffuse radiation - early Universe 


1 INTRODUCTION 


CMB temperature, T 7 (e.g. Furlanetto et al. 2006): 


Radiation emitted from the first stars and galaxies marked 
the end of the cosmic dark ages. The UV ionizing photons 
from these first sources escape into the intergalactic medium 
(IGM), reionizing the pervasive neutral hydrogen fog, bring¬ 
ing about the final major baryonic phase change of the Uni¬ 
verse. This epoch of reionization (EoR) is rich in information 
on the formation and evolution of structures in the early 
Universe, including the nature of the first stars and galax¬ 
ies and their impact on the IGM (see e.g. Barkana & Loeb 
2007; Loeb & Furlanetto 2013; Zaroubi 2013). 

The most promising observational tool for probing the 
EoR physics is the redshifted 21-cm spin-flip transition of 
neutral hydrogen(see e.g. Gnedin & Ostriker 1997; Madau 
et al. 1997; Shaver et al. 1999; Tozzi et al. 2000; Gnedin & 
Shaver 2004; Furlanetto et al. 2006; Morales & Wyithe 2010; 
Pritchard & Loeb 2012). The signal is expressed as the offset 
of the 21-cm brightness temperature, 5Tb(v), relative to the 


5Tb (n) 
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where xm is the neutral fraction, Ts is the gas spin tempera¬ 
ture, <5 n i(cc,2) is the evolved (Eularian) overdensity, H(z) is 
the Hubble parameter, du r /dr is the gradient of the line-of- 
sight component of the velocity and all quantities are eval¬ 
uated at redshift z = v o/V — 1, where uo is the 21-cm fre¬ 
quency. Thus, the 21-cm signal probes both the ionization 
and thermal state of the IGM, and indirectly also the galax¬ 
ies which regulate them. Moreover, as it is a line transition, 
it promises a 3D view of the early Universe. 

Tapping into this physical bounty requires sensitive ra¬ 
dio interferometers. First generation experiments such as the 
Low Frequency Array (LOFAR; van Haarlem et al. 2013; 
Yatawatta et al. 2013) 1 , the Murchison Wide Field Array 
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1 http://www.lofar.org/ 
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(MWA; Tingay et al. 2013) 2 and the Precision Array for 
Probing the Epoch of Reionization (PAPER; Parsons et al. 
2010) 3 seek the statistical characterization of the 21-cm sig¬ 
nal through a power spectrum (PS) measurement. Though 
no detection has been forthcoming thus far, steady progress 
is being made. Presently, the tightest upper limits come from 
the 64-dipole PAPER array, which constrained the spheri¬ 
cally averaged PS over 0.1 < k < 0.4 h/ Mpc to be less than 
< 500 mK 2 . Tomographic maps will have to wait for second- 
generation instruments such as the Square Kilometre Array 
(SKA; Koopmans et al. 2015) 4 and the Hydrogen Epoch of 
Reionization Array (HERA; Beardsley et al. 2015) 5 . 

Even these preliminary observations begin to constrain 
the state of the early Universe. From equation (1), we see 
that, once the IGM spin temperature is much hotter than the 
CMB (Ts T 7 ), the maximum achievable contrast in the 
21-cm signal is of order 10 mK, sourced by order unity fluc¬ 
tuations in the neutral fraction during patchy reionization. 
However, prior to IGM heating when the signal is still in 
absorption against the CMB, the achievable dynamic range 
is much larger. For example, the adiabatically cooling IGM 
can reach values of (5Tb ~ —200mK. 6 7 Although the spatial 
fluctuations in the IGM temperature source a strong 21-cm 
PS signal (~ 100s of mK 2 ; e.g. Pritchard & Furlanetto 2007; 
Pacucci et al. 2014), an even stronger signal would arise 
during so-called cold reionization (Mesinger et al. 2013). If 
reionization proceeded in a very cold IGM, the resulting 
contrast between the cold neutral regions (with <5Tb ~ —100 
to —200 mK), and the ionized regions (with 5Tb ~ 0 mK) 
could drive the large-scale 21-cm power to values in excess 
of ~ 1000s of mK 2 (e.g. Parsons et al. 2014). Such models 
can already be ruled out by current data. 

Parsons et al. (2014) were the first to investigate con¬ 
straints on ‘cold reionization’ from upper limits on the 
z = 7.7 21-cm PS using the 32 dipole PAPER array. Ali 
et al. (2015), using tighter upper limits on the 21-cm PS at 
2 = 8.4 from an extended 64 element PAPER array, subse¬ 
quently obtained improved limits of Ts > 4 K. Both works 
only used a simple analytic expression for the 21-cm PS dur¬ 
ing the EoR. An improved analysis was performed in Pober 
et al. (2015), making use of the semi numerical code 21CM- 
FAST' (Mesinger & Furlanetto 2007; Mesinger et al. 2011) 
to simulate both the EoR and IGM heating by X-rays. These 
authors varied the X-ray efficiency of early galaxies, in or¬ 
der to coarsely sample the IGM spin temperature, obtaining 
lower limits of Ts ^ 5 K for neutral fractions between 10 and 
85 per cent. 

Importantly, in each of these studies, only a single EoR 
model was used. The large-scale morphology (distribution of 
cosmic H ii patches) depends on the properties of the sources 
and sinks of ionizing photons. These can produce large vari¬ 

2 http://www.mwatelescope.org/ 

3 http://eor.berkeley.edu/ 

4 https://www.skatelescope.org 

5 http://reionization.org 

6 This signal is achieved if the spin temperature is efficiently 
coupled to the gas kinetic temperature, through the Won I huysen- 
Field mechanism (Wouthuysen 1952; Field 1958). Simple scaling 
relations (e.g. Furlanetto 2006; McQuinn 2012) show that this 
condition should be satisfied well before the IGM heating epoch. 

7 http://homepage.sns.it/mesinger/Sim 


ations in both the shape and amplitude of the 21-cm PS 
during the EoR, even at a fixed redshift and mean neutral 
fraction (e.g. fig. 2 of Greig & Mesinger 2015). 

In order to explore the impact of the EoR morphology 
on the constraints from PAPER-64, we utilize the recently 
developed Monte Carlo Markov Chain (MCMC)-based EoR 
analysis tool 21CMMC 8 (Greig & Mesinger 2015). Using this 
theoretical framework, we recover improved, robust lower 
limits on the IGM temperature by marginalising over EoR 
models. We then strengthen these limits, by including reion¬ 
ization priors from observations of the dark pixel statistics of 
high-z quasars (McGreer et al. 2015) and the electron scat¬ 
tering optical depth, r e (Planck Collaboration XIII 2015). 

The remainder of this paper is organized as follows. In 
Section 2, we outline our analysis methodology, while in Sec¬ 
tion 3, we discuss our improved constraints on the IGM tem¬ 
perature. Finally, in Section 4, we finish with our closing re¬ 
marks. Throughout this work, we adopt the standard set of 
ACDM cosmological parameters: (O m , Q,a, fib, n, as, Ho) = 
(0.27, 0.73, 0.046, 0.96, 0.82, 70 kms _1 Mpc _1 ), measured 
from WMAP (Bennett et al. 2013) which are consistent with 
the latest results from Planck (Planck Collaboration XVI 
2014). 


2 METHODOLOGY 

In Greig & Mesinger (2015), we developed 21CMMC, an 
MCMC-based analysis tool enabling the exploration of the 
EoR astrophysical parameter space. To simulate the reion¬ 
ization epoch, we assume a relatively popular, three param¬ 
eter model which can accommodate a large set of physically 
motivated EoR morphologies. In this section, we summarize 
the EoR model sampled within 21CMMC and the modifi¬ 
cations required to model ‘cold reionization’, deferring the 
reader to Greig & Mesinger (2015) for more technical dis¬ 
cussions. 


2.1 Modelling reionization within 21CMMC 

For a given EoR parameter set, we simulate the 21-cm sig¬ 
nal using the semi numerical simulation code 21CMFAST. 
21CMFAST employs approximate but efficient methods for 
modelling the 21-cm signal, and is accurate when com¬ 
pared to computationally expensive radiative transfer simu¬ 
lations on scales relevant to 21-cm interferometry, ^ 1 Mpc 
(Mesinger et al. 2011; Zahn et al. 2011). The speed and ef¬ 
ficiency of 21CMFAST makes it well suited for MCMC sam¬ 
pling. Below we outline the basic components of the EoR 
simulation; readers interested in more details are encour¬ 
aged to read Mesinger & Furlanetto (2007) and Mesinger 
et al. (2011). 

21CMFAST produces 3D realizations of the IGM den¬ 
sity, velocity, source and ionization fields. A cubic volume 
of the linear density field is evolved using the Zel’dovich ap¬ 
proximation (Zel’dovich 1970). The ionization fields are esti¬ 
mated using the excursion-set formalism outlined in Furlan¬ 
etto et al. (2004), but modified to operate on the evolved 
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density field (Mesinger et al. 2011). Following this prescrip¬ 
tion, the time-integrated number of ionizing photons are 
compared to the number of baryons within spherical regions 
of decreasing radius, R. A cell within the simulation volume 
is then classified as fully ionized if, 

Cfcon(x,z, R, M min ) ^ 1, (2) 

where £ is the ionization efficiency which describes the con¬ 
version of mass into ionizing photons (see Section 2.1.1) and 
fcoii(x, z, R, M m i n ) is the fraction of collapsed matter within 
a spherical radius R residing within haloes larger than M m i n 
(Press & Schechter 1974; Bond et al. 1991; Lacey & Cole 
1993; Shetli & Tormen 1999). Any cells not fully ionized 
(partial ionizations) are smoothed by the minimum smooth¬ 
ing scale of the cell, R ce n, and their ionization fraction from 
internal sources is set to Cf co n{x, z, R ce n, Mmin). 

In this work, we adopt a popular three parameter model 
to characterize the EoR: (i) the ionizing efficiency of high -2 
galaxies; (ii) the mean free path of ionizing photons; (iii) the 
minimum virial temperature hosting star-forming galaxies. 
These EoR parameters are somewhat simplistic as they in 
effect average over redshift and/or halo mass dependences. 
However, this simple model suffices to describe a broad range 
of EoR morphologies, while at the same time providing a 
straightforward physical interpretation. In addition to these 
EoR parameters, here we include an additional free param¬ 
eter: the mean IGM spin temperature. To help gain some 
intuition, below we summarize these four parameters, high¬ 
lighting how they affect the 21-cm signal. 


2.1.1 Ionizing efficiency, £ 


The ionizing efficiency of high -2 galaxies (equation 2) can 
be expressed as 


C = 30 




N, \ 
4400 ) 


1.5 


1 + n r . 


(3) 


where, / esc is the fraction of ionizing photons escaping into 
the IGM, /* is the fraction of galactic gas in stars, AQ is the 
number of ionizing photons produced per baryon in stars and 
nrec is the average number of recombinations per baryon in 
the IGM. 

The parameter £ mainly serves to speed up/slow down 
reionization. Within this work, we take a flat prior over the 
range £ £ [5,200], which results in a range of reionization 
histories which are in broad agreement with current EoR 
constraints (Greig & Mesinger, in prep). 


where p is the mean molecular weight, 17]/ = I7 m (l + 
2 ) 3 /[fi m (l + z) A + 17a], and A c = 187t 2 + 82 d — 39 d 2 where 
d = f7m — 1. The value T™ r m ~ 10 4 K corresponds to the 
minimum temperature for efficient atomic cooling; however, 
efficient star formation likely requires more massive haloes, 
which are able to better retain and reincorporate supernovae 
driven outflows (e.g. Springel & Hernquist 2003). 

The parameter T™ r m affects (i) when (at what redshift) 
reionization occurs, and (ii) the bias of the galaxies which are 
responsible. A higher value of T™ r m means that reionization 
happened later, with more large-scale ionization structure 
(at a fixed value of the mean neutral fraction; e.g. McQuinn 
et al. 2007). Here we assume a flat prior over the log of 
the virial temperature, within the range T™ r m £ [10 4 ,5 x 
10 s ] K. The lower limit corresponds to the atomic cooling 
threshold and the upper limit is roughly consistent with the 
host haloes of observed Lyman break galaxies at 2 ~ 6-8 
(e.g. Kuhlen & Faucher-Giguere 2012; Barone-Nugent et al. 
2014). 


2.1.3 Mean free path of ionizing photons within ionized 
regions, R m f P 

The physical scales on which escaping ionizing photons 
are able to penetrate into the surrounding IGM depends 
strongly on both the number and properties of the absorp¬ 
tion systems (such as Lyman limit and more diffuse sys¬ 
tems). Typically below the resolution limits of EoR simula¬ 
tions, these systems behave as photon sinks, resulting in a 
maximum physical scale on which H ii regions are capable of 
growing around the ionizing galaxies. The impact of photons 
sinks can be crudely interpreted as a maximum horizon for 
the ionizing photons, which can correspond to the maximum 
filtering scale in excursion-set EoR models. We denote this 
scale as 7? m f p , noting that it depends on the time-integrated 
value of the mean free path through the ionized IGM during 
the EoR (Sobacchi & Mesinger 2014). 

The parameter i? m f P most strongly impacts the large- 
scale ionization structure, with the large-scale 21-cm power 
dropping with decreasing i? m f p . Here, we adopt a flat prior 
over the range R m i P £ [5,40] cMpc, motivated by sub¬ 
grid models of inhomogeneous recombinations (Sobacchi & 
Mesinger 2014), as well analytic estimates (Furlanetto & 
Oh 2005) and hydrodynamical simulations of the IGM (Mc¬ 
Quinn et al. 2011; Emberson et al. 2013). 


2.1.2 Minimum virial temperature of star-forming haloes, 

ry-imin 

vir 


The minimum threshold for a halo hosting a star-forming 
galaxy, regulating processes important for star formation 
such as gas accretion, cooling and retainment of supernovae 
outflows, can be defined in terms of its virial temperature, 
T™ r m , which is related to its halo mass via, (e.g. Barkana & 
Loeb 2001) 


M n: 


= 10 s h 


V0.6 ) 


-3/2 / f7 


- 1/2 


T vi 


1.98 x 10 4 K 


f7; 

3/2 


187t 2 / 

1+2 

10 


-3/2 


M 0 , (4) 


2.1.4 Mean IGM spin temperature, Ts 

In order to recover limits on the ‘cold reionization’ model, 
in this work, we adopt a global value for the IGM spin tem¬ 
perature of the neutral cosmic gas during the EoR. This is 
an oversimplification, as we would expect spatial fluctua¬ 
tions in Ts sourced by the clustering of X-ray sources and 
the corresponding mean free path of the X-rays. However, 
in the relevant regime when the signal is the strongest (ad¬ 
vanced stages of reionization with weak X-ray heating), a 
uniform Ts is a good approximation for most of the IGM as 
the hottest parts of the IGM which normally source large 
spatial fluctuations in Ts have already been ionized. The 
remaining neutral gas is distant from galaxies and its heat¬ 
ing is governed by the harder X-rays with longer mean free 


© 0000 RAS, MNRAS 000, 000-000 



4 B. Greig et al. 


paths. This picture was quantitatively confirmed by Pober 
et al. (2015), who note that ignoring the spatial fluctuations 
in the spin temperature of the neutral IGM only impacts 
the large-scale 21 -cm power during ‘cold reionization’ at the 
level of ~ 10 per cent (compared with the factor of 10 de¬ 
pendence of the large-scale power on the EoR morphology; 
e.g. fig. 2 in Greig & Mesinger 2015). 

Within this work we explore the range Ts € [1.9, 30] K. 
This range encompasses the low IGM temperatures required 
to source a 21 -cm signal strong enough to be ruled out by the 
upper limits in Ali et al. (2015). For reference, the IGM tem¬ 
perature in the absence of heating by astrophysical struc¬ 
tures at 2 = 8.4 corresponds to 1.9 K (RECFAST; Seager 
et al. 1999, 2000), whereas the CMB temperature at that 
redshift is 25.6 K. 

2.2 Computing the likelihood with 21CMMC 

The principal data used for this analysis are the PS measure¬ 
ments from PAPER-64 presented in Ali et al. (2015). Here, 
we summarize the significant properties of their analysis; re¬ 
ferring the reader to their work for detailed discussions of 
the observations and data reduction. Ali et al. (2015) anal¬ 
ysed the complete observational campaign conducted with 
a 64-element PAPER array in South Africa over the span of 
135 days from 2012 November to 2013 March. Their analysis 
focused on a 10 MHz frequency band centred on 151.5 MHz, 
corresponding to a 21-cm redshift of 2 = 8.4. The end re¬ 
sult is a quoted upper limit on the 21-cm PS of (22.4 mK ) 2 
over the range of 0.15 < k < 0.5 h/ Mpc, nearly a factor of 
four lower (in mK 2 ) than the previous best upper limit of 
Parsons et al. (2014). Approximately half of this increased 
sensitivity comes from a doubling of the number of antennas 
(64 elements, as opposed to the 32 used in the Parsons et al. 
2014 analysis). The other half stems from an updated anal¬ 
ysis method, with the principal three advances coming from 
improved redundant calibration using the methods of Zheng 
et al. (2014), fringe rate filtering (Parsons et al. 2015), and 
optimal quadratic estimators (Liu & Tegmark 2011; Trott 
et al. 2012; Dillon et al. 2013; Liu et al. 2014a,b). 

The final PS measurements of Ali et al. (2015) corre¬ 
spond to nine data points spanning the range 0.1 < k < 
0.5 h/ Mpc, and their corresponding 2a errors were obtained 
from a boot-strapping analysis. Within this work, we con¬ 
sider only data points with 2 a errors that are consistent with 
a null detection, as systematic errors remain in the data and 
will lead to significant excesses of power. Indeed within Ali 
et al. (2015), the ‘detections’ at low k are reported as likely 
being a result of foregrounds, whereas for the other ‘detec¬ 
tions’ the causes are less certain. As a result, we are left 
with a total of four useable data points, which we can set as 
upper limits on the shape and amplitude of the 21-cm PS. 

In the introductory paper of Greig & Mesinger (2015), 
we performed a maximum likelihood sampling of the EoR 
parameter space using a y 2 statistic. To incorporate only 
the upper limits on the 21-cm PS, we retain the y 2 statistic, 
however, we modify the likelihood computation. If the am¬ 
plitude of the model PS falls below the observed value, we 
set our x 2 to zero, resulting in equal likelihoods for the cor¬ 
responding EoR parameters. If instead the model amplitude 
is above the observation, we compute the standard y 2 , using 
the bootstrapped errors from Ah et al. (2015). Such an ap¬ 


proach ensures we appropriately down weight EoR models 
producing amplitudes brighter than the PAPER-64 limits, 
while accepting all models fainter 9 . 

Finally, we remove the conservative 25 per cent mod¬ 
elling uncertainty adopted in Greig & Mesinger (2015) when 
simulating the PS from 21 CMFAST. This source of error 
is negligible when compared with the high upper limits on 
the 21-cm PS. Furthermore, we remove our conservative k- 
mode cut at k = 0.15 Mpc -1 , which was included to remove 
Fourier modes contaminated by foregrounds, since this was 
already accounted for in the PAPER-64 limits. 


2.3 Including observational priors on the EoR 

In an effort to further improve the constraints on the IGM 
temperature, we consider two observational priors. First, we 
consider the strongest available constraints on the tail of the 
reionization epoch through measurements of the dark-pixel 
fraction of high -2 quasars (McGreer et al. 2015). This ap¬ 
proach provides model independent constraints on the IGM 
neutral fraction at 2 = 5.9 of xm ^ 0.06 + 0.05 (la), indica¬ 
tive of a completely (or almost nearly) reionized Universe 
by 2 « 6 . 

Secondly, we consider the measurement of the electron 
scattering optical depth to the CMB, r e = 0.066±0.016 (la) 
(Planck Collaboration XIII 2015). This corresponds to an 
instantaneous reionization redshift of z le = 8 . 8 j 4 ' 4 , close 
to the observational redshift of the PAPER-64 21-cm PS 
constraints (2 = 8.4). 

These priors are incorporated into 21CMMC by com¬ 
puting the 21-cm PS at five different epochs, 2 = 

6,7,8.4,10 and 11 for each EoR model parameter set in 
the MCMC, whereas our y 2 estimate is only performed at 

2 = 8.4. We confirm that this sampling of the reionization 
history ensures we obtain accurate estimates for r e . 

3 RESULTS 

3.1 Constraints only from PAPER—64 

In the top-left panel of Fig. 1, we report the joint 2D likeli¬ 
hood constraints of the IGM temperature and the IGM neu¬ 
tral fraction recovered from 21CMMC. To model the prob¬ 
ability surface, we assume a normalized Gaussian distribu¬ 
tion whereby each point in the Ts-*hi plane is assigned the 
probability, P(a;Hi,Ts) oc exp(—|y 2 ), where the y 2 value is 
determined by 21CMMC. In order to obtain this parameter 
surface, we marginalise over our three morphological EoR 
parameters (£, R m f p and T™ r m ) assuming flat (uniform) pri¬ 
ors for each across their allowed astrophysical values (see 


9 This approach differs to that adopted by Pober et al. (2015). 
These authors use all data points as an upper limit, not just 
those consistent with zero. For each fc-mode, they then compute 
the probability of obtaining their model 21-cm PS relative to the 
observation, and construct the likelihood by determining the joint 
probability of obtaining all data points for the model 21-cm PS. 
In doing so, these authors allow slightly higher amplitude 21-cm 
PS models than our approach: however, the relative differences 
between the two approaches should be minimal. 
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Neutral Fraction, xh i Neutral Fraction, cch i 

Figure 1. Joint 2D likelihood surfaces for the IGM spin temperature, Tg and the IGM neutral fraction, Sjji> marginalised over our 
three morphological EoR parameters (£, ft m f p and T™ r m ) assuming flat (uniform) priors. Shaded regions correspond to P(xu \, Tg) J 0.6 
(approximating ltr), while dashed lines correspond to P p 0.15 (approximating 2cr). In the top-left panel, we provide the lower limits 
obtained from only the PAPER-64 constraints on the 21-cm PS. Top right, we include a prior on the IGM neutral fraction at z = 5.9 
(,T|.|i ± 0.06 + 0.05, McGreer et al. 2015) to the PAPER-64 data. Bottom left, the inclusion of a prior on the electron scattering optical 
depth, T e (r e = 0.066±0.016, Planck Collaboration XIII 2015). Bottom right, our tightest constraints on the IGM temperature, including 
both the ini[ z = 5.9] and r e constraints. In all panels, the grey shaded region represents the lower limit on the IGM temperature in the 
absence of heating by astrophysical structures, Tg = 1.9 K. 


Section 2.1). In all panels of Fig. 1, the shaded regions cor¬ 
respond to likelihoods of P{xm,Ts) P 0.6; we consider the 
regions in white, corresponding to P(xm,Ts) ^ 0.6, to be 
disfavoured by the data 10 

From the top-right panel, we see that the current 
PAPER-64 limits constrain the IGM temperature to be 
Tg > 5 K at greater than 60 per cent confidence, for an 

10 Since the observation only consists of an upper limit on the 
21-cm PS, the values of Tg —> oo are allowed with relatively equal 
likelihood (see equation 1). Therefore, constructing standard con¬ 
fidence limits on Tg (defined to enclose a fixed area of the prob¬ 
ability density function) is sensitive to the chosen range of the 
Tg priors. We instead choose to present our results as contours of 
equal likelihood, which are well defined, although the exact choice 
of the threshold value is arbitrary. The fiducial choice of P p 0.6 
is loosely motivated by the probability of a normalized Gaussian 
distribution at 1 a. For illustrative purposes, in all panels in Fig 1, 
we additionally highlight an adopted threshold of P P 0.15 cor¬ 
responding to the probability of a normalized Gaussian at 2a. 


IGM neutral fraction between 30 and 65 per cent. These re¬ 
duce to Ts > 3 K for a broader range in the IGM neutral 
fraction of xm ->10 per cent. 

Qualitatively, the shape of the allowed Tg-aSm param¬ 
eter space is almost identical to that presented in Pober 
et al. (2015). As discussed in that paper, this shape arises 
from the fact that the spin temperature term in equation 
(1) can roughly act as a multiplicative factor of the 21-cm 
PS during the EoR under the usual Tg T 7 assumption. 
This PS peaks during the midpoint of reionization, and in 
that regime the IGM does not need to be as cold to exceed 
the Ali et al. (2015) upper limits. 

Quantitatively however, our lower limits on the IGM 
spin temperature are reduced by at least a factor of two 
compared to those in Pober et al. (2015). It is difficult to 
pin down all causes for this discrepancy, owing to the afore¬ 
mentioned different approaches in constructing the likeli¬ 
hood. Nevertheless, a relaxed lower limit is consistent with 
our naive expectations. The single EoR model used in Pober 
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et al. (2015) results in sizeable ionization structure on large 
scales. On the other hand, our EoR parameter space addi¬ 
tionally includes models with substantially less large-scale 
ionization structure (compared at a fixed neutral fraction; 
see fig. 2 in Greig & Mesinger 2015), mostly sourced by 
smaller mean free paths of ionizing photons (e.g. Sobacchi 
& Mesinger 2014). Therefore, by marginalizing over all EoR 
models within our MCMC framework, we would expect a 
loosening of the lower limits owing to the broadening of the 
recovered probability distribution for the IGM temperature. 

3.2 Constraints including additional EoR priors 

In the second panel of Fig. 1, we include constraints on the 
end of reionization at 2 « 6 from McGreer et al. (2015). We 
find this drastically reduces the allowed parameter space, 
disfavouring models with xm > 0.75. This behaviour is 
straightforward to interpret. To produce a ~ 10 per cent 
neutral IGM by 2 « 6 it would be very difficult to have a 
large aim at 2 = 8.4, requiring unphysically-rapid reioniza¬ 
tion histories. 

Next, we consider the inclusion of the r e prior from 
Planck , shown in the bottom-left panel of Fig. 1. The Planck 
measurement favours a midpoint of reionization close to the 
same redshift as the PAPER-64 measurement (albeit with 
broad error bars). This disfavours models with either too 
low or too high a neutral fraction, which are precisely those 
in which the signal is least sensitive to Ts. Thus our lower 
limits on Ts improve with respect to the top panels. 

Finally, in the bottom right panel, we provide our lower 
limits combining both observational priors. Owing to the 
complimentary nature of the two priors, as discussed above, 
we now recover our most stringent lower limits on the IGM 
temperature of Ts > 6. Even with a lower limits of Ts « 6 K, 
we can begin to place constraints on reionization scenar¬ 
ios such as ‘cold reionization’. As this temperature is above 
that of the IGM in the absence of heating by astrophysi- 
cal structures (1.9 K), this implies that the IGM must have 
undergone some level of heating, ruling out a truly ‘cold’ 
(unheated) IGM. However, the relatively small amount of 
heating required does not yet begin to strongly constrain 
the physical processes that could be responsible (e.g. Pacucci 
et al. 2014). 

Nevertheless, this approach highlights the utility of cur¬ 
rent 21-cm experiments. As more data are acquired, with 
better characterized systematics, these upper limits on the 
21-cm PS amplitude should reduce, tightening the lower lim¬ 
its on the IGM temperature. Alternately, should a high 21- 
cm PS actually be detected, the implied lack of X-ray heat¬ 
ing would have significant consequences for the star forma¬ 
tion physics of the first galaxies. 


4 CONCLUSION 

One of the major astrophysical goals for the forthcoming 
decade is the statistical measurement of the 21-cm PS from 
the EoR. In probing the phase change from a neutral to ion¬ 
ized IGM, the EoR will yield rich information on the forma¬ 
tion, growth and evolution of the first stars and galaxies and 
their influence on the IGM. For this reason several, dedicated 
first-generation radio instruments have been constructed, to 


be followed by a wave of second-generation experiments with 
significantly larger collecting areas and sensitivity. 

Although these experiments have yet to provide a de¬ 
tection, even current upper limits on the amplitude of the 
21-cm PS can begin to constrain the physics of the EoR. In 
the absence of X-ray heating (or any other heating mecha¬ 
nism), the neutral IGM should remain ‘cold’, well below the 
CMB temperature. If reionization proceeds in such a cold 
IGM, the resulting contrast between cosmic ionized patches 
(with 5Tb ~ 0 mK) and the cold, neutral patches (with 
5Tb ~ —100 mK) would drive the 21-cm signal to values 
in excess of current upper limits. Recently, both Ali et al. 
(2015) and Pober et al. (2015) used this fact to place limits 
on the IGM temperature during reionization, based on the 
21-cm PS upper limits with PAPER-64. 

In this work, we improve upon these existing analyses 
in several ways. Using a modified version of the MCMC- 
based EoR analysis tool 21CMMC (Greig & Mesinger 2015) 
we better sample the allowed astrophysical parameter space, 
improving on the fixed grid approach of Pober et al. (2015) 
and the analytic ‘toy’ model of Ali et al. (2015). This is 
achieved by allowing the mean IGM spin temperature, Ts 
to be a free parameter rather than varying the X-ray effi¬ 
ciency to recover a mean Ts as in Pober et al. (2015). Most 
importantly, by employing 21CMMC, we marginalize over a 
broad range of EoR models and morphologies. Both previous 
analyses assumed a single EoR model, thereby neglecting the 
fact that the EoR morphology can impact the amplitude of 
the 21-cm PS by up to an order of magnitude at a fixed 
neutral fraction (e.g. Greig & Mesinger 2015). 

From the PAPER-64 upper limits alone, we can recover 
lower limits on the IGM temperature at 2 = 8.4 of Ts > 5 K 
at greater than 60 per cent confidence for an IGM neutral 
fraction between 30 < Shi < 65 per cent, which reduces 
to T S > 3 K when xm > 10 per cent. We further tighten 
these constraints by including EoR priors. Folding-in priors 
on the tail end of the reionization epoch (aim ^ 0.06 + 0.05 
at 2 = 5.9; McGreer et al. 2015), as well as the integral 
constraint of r e = 0.066 ± 0.016 (Planck Collaboration XIII 
2015), we obtain lower limits of Ts > 6 K at greater than 60 
per cent confidence. While not completely ruling out ‘cold 
reionization’ models, our results confirm that at least some 
level of prior IGM heating occurred. As both the quality 
and volume of available 21-cm data continues to increase, 
our framework will be able to provide unique insights into 
the high-energy processes inside the first galaxies. 
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